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Abstract 

An off-line version of the Wall Interference Correc- 
tion System (WICS) has been implemented foi the 
NASA Langley National Transonic Facility. The cor- 
rection capability is currently restricted to corrections 
for solid wall interference in tin* model pitch plain 1 for 
Mach numbers less than 0.45 due to a limitation in 
tunnel calibration data. A study to assess output sen- 
sitivity to the aerodynamic parameters of Reynolds 
number and Mach number was conducted on this code 
to further ensure quality during the correction process. 
In addition, this paper includes an investigation into 
possible correction due to a semispan test technique 
using a non metric standoff and an improvement, to 
the standard data rejection algorithm. 


Nomenclature 


Cp 

Nondimensional drag coefficient 

Ci. 

Noiidimensional lift coefficient 

M 

Mach number 

P 

Pressure 

Re 

Reynolds number per unit length 

U v 

Axial velocity determined by calibration 

Etc 

Free-stream velocity 

i/M M Wi 

Velocity perturbation: .r, y, z components 

!/, ~ 

Cartesian body axes 

n 

Angle of attack 


Angle of sideslip 

A 

Correction 
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^ Velocity potential 

Subscripts 
row Corrected 

F Free air 

j incremental correction 

so scmispa.ii standoff 

T In tunnel 

unr Uncorrected 

oc Freestr earn 

Introduction 

The subject of wall interference has broadened im- 
mensely over the past century. Classical correction 
methods for this effect incorporate a mathematk al 
boundary condition for the wind tunnel walls. Non- 
dimensional parameters have been tabulated and a 
number of formulas to correct spatially- varying wall- 
interference for aircraft-like models have been docu- 
mented in texts such as AGARDograph 109. 1 Moie 
often facilities arc' testing in regions that exceed the 
design specifications of the facility. Increased model 
size, for instance, may introduce larger amounts of 
blockage and lift interference than the classical meth- 
ods can properly handle. Modern, more 1 capable 
interference-correction methods, including boundary 
pressure methods were recently published in AGAR- 
Dograph 336. 2 The advantage of the boundary pres- 
sure methods over classical methods is their ability 
to respond to the actual conditions of the flow in the 
tunnel test section; additionally, boundary pressure 1 
methods are more complex than classical methods but 
can be implemented for online post-point /post-scan 
applications. 

Customers arc 1 also placing more 1 stringent accuiac\ 
requirements on the acquired data. To reduce the 
uncertainties introeiuced in tunnel data due to the 
presence e)f the walls, thereby enhancing data quality. 
NASA Langley Research Center has implemented the 
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Wall Interference Correction System (WIGS) for solid 
wall testing in the National Transonic Facility (NTF). :i 
The code was originally developed hy Ulbrich at al x u 
at NASA Arnos Research Center for the 12-Foot Pres- 
sure' Wind Tunnel. The W IC S code' uses a modified 
Hackett (Wall Signature) Method to eletermine wall 
interference corrections. This method was selected be- 
cause it is rapiel and robust, and it has minimal impact 
on the facility in terms of instrumentation when com- 
pared to other boundary pressure 1 methods. A brief 
description of the theory behind the WTCS method 
and the implementation at the NTF is presented. 

The purpose of this paper is to determine cock 1 out- 
put sensitivities to aerodynamic parameters, support 
systems, and other input, parameters for both semispan 
and fullspan models. These cases will be summarized 
for application during operational use 1 of the WTCS 
code. 

Correction Accuracy Requirements 

Before a sensitivity study can he conducted on the 
WTCS code, it is necessary to understand the accuracy 
requirements placed on wall interference corrections in 
general. According to Steinlc and Stanewsky, 10 the 
requirements for transonic cruise are dominated by 
a one count (Co = 0.0001) accuracy for drag for a 
transport-type aircraft. Since the WTCS code is not 
used exclusively on transport-type 1 aircraft in cruise, a 
more general requirement is needed. Newman and Ev- 
erhart 2 report the accuracy requirements listed below 
in Table 1. 


Type of Test 

Incremental 

Absolute 

High Lift 

0.2% O l A' C u 

0.4% C L kC D 

Transonic 

1/2 count Cd 

1 count Co 


Table 1 Generalized Industry Correction Accu- 
racy Requirements 


Description of the NASA LaRC 
National Transonic Facility 

The NTF 11 is a fan-driven, closed-circuit, 
continuous- flow, pressurized wind tunnel, which 
is capable of testing at cryogenic conditions. The 
test gas is dry air for elevated temperature operation 
and nitrogen for reduced temperature operation. 
The settling chamber contains four anti-turbulence 
screens. A 15:1 contraction ratio entrance cone leads 
into an 8.2 feet square cross sectional test section with 
six inch triangular corner fillets which extends 25 feet 
in length then opens into a high speed diffuser. The 
operational envelope of the NTF encompasses a large 
range of test conditions. The facility can sustain a 
continuous airspeed from 0.1 to 1.2 in Mach number. 
Total pressure capabilities of the facility range from 
15 to 130 psi. The tunnel can operate at temperatures 
ranging from 150° F down to -320°F. These large 


ranges of conditions allow Reynolds number testing 
from 3 to 120 million per foot, NTF has the capa- 
bility to independently vary Mach number, Reynolds 
number, and dynamic pressure. 

Both fullspan and semispan model mounting sys- 
tems are available in the NTF. Fullspan models are 
supported by an aft mounted sting. This sting is 
attached to a vertically oriented arc sector used to 
change model pitch attitude over a range from -11° 
to 19°. The roll mechanism has a range of -90° to 
180". Sideslip angles are achieved by using combined 
roll and pitch angles. A sidewall mounting system is 
used for semispan models. The 1 angle of attack range 
is ±180°. The center of rotation for semispan and 
fullspan models is tunnel station 13 feet. 

Currently, the NTF has 459 operational wall pres- 
sure orifices of which the WTCS code uses 360. Looking 
downstream. Figure 1 shows the? cross sectional dia- 
gram of the rows of pressure ports; it also serves as 
a map between the NTF row numbering scheme and 
that used for WTCS. Figure 2 unwraps the tunnel walls 
to show the current wall orifice layout. The filled port 
symbols are an example subset of port selections used 
for WTCS. 2 

The Wall Interference Correction 
System (WICS) 

The following presents a brief overview of Ulbrich \s 
extension of the Hackett wall signature method, in- 
cluding theory and the current implementation of 
W ICS in the NTF. A more detailed description is 
given by Ulbrich et al and Iyer et al A The 
method applies a measured pressure boundary condi- 
tion which is the fared difference 1 between the model 



Fig. 1 Cross-Sectional Pressure Measurement 
Row Setup of the NTF 
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Fig. 2 Wall Orifice Layout for NTF 


installed condition and the empty tunnel. The test 
article is comprised of singularities: a source-sink com- 
bination is used to represent the fuselage: sources art' 
used to represent the viscous separation in the wake; 
and lint* doublets, typically placed on a lifting surface* 
quarter-chord, art* used to simulate the effect of lift. 
Additionally, a powered simulation can he represented 
by proper placement of a sink. Once the strengths 
of the line doublets and propulsion simulation sinks 
arc* known, the* signature of these singularities is sub- 
tracted from the t. arc'd wall signature*. The remaining 
signature* is used to determine the strengths of the solid 
and wake blockage singularities. Hackett et ai l2 ] * 
originally used a “local” least squares fitting procedure 
of the wall pressure data, which was labor intensive 
because each pressure* measurement used in the fit of 
the* wall signature had to be individually selected for 
each data point. In an effort to improve* the process. 
Ulbrich 9 introduced a “global” least squares fitting 
procedure* that uses all available pressure ports. The* 
use* of panel code's to precompute normalized solutions 
for use* in matching the* signature was also introduced. 
Ulbrich further modified the Hackett method by us- 
ing balance measured force's and moments to calculate* 
the* strength of the doublet singularities. The method 
use*s tareel data to assess the wall induced effects on 
the test article. The* t.ared data are* obtained by sub- 
tracting empty tunnel wall pressure port calibrations 
from the* moeleTinstalle*d pressure data to remove first 
order effects of the* empty tunnel hemndary layer and 
buoyancy. This taring of the measured data assumes 
that the* additional second e>rde*r change in the bound- 


ary laye*r displacement thickne'ss due te) the* model 
in the tunnel is negligible*. Normalized perturbation 
velocities are calculated from the* subsonic potential 
equation using the* method of images. fe>r each type* 
of singularity in discrete locations in the* tunnel. By 
superimposing these* model singularities with the'ir e'en- 
responding calculated strengths, interference' velocities 
are computed. 

The* blejckage inte*rfe*re*nce factor c is defined by 


Hi ^ ur - ur 


(i) 


wher i'Ui/U,. = dipr/dx nm\(pj « Incidence 

and sideslip corrections are de*fined by 


Oti 

Pi 


U< tt'j ^ U t wr - u'r 
U ( r, ^ l ( vt ~ v i' 

~ u ~ u r “ 


( 2 ) 

( 3 ) 


where* Uy^/U ( is elefine*d by 1 + s. It is important 
to note* that, the* intograteel e*ffe*ct of tlie* e*mpty tun- 
nel boundary layer growth and buoyancy due to im- 
proper wall divergence arc* alre*ady iiiclude*d in l , . The* 
“empty-tunnel” calibration provides the* perturbation 
velocity elenoted by a*. which e*an be used to correct 
the measured moele*Tin-tunnel perturbation velocity. 
ulj . by 

ur Ur - us /n 

= tv !l: 

The ratio wj /U, is used to determine* the* singularity 
strengths. 
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support system kinematics must, also be considered. 
Semispaii and fullspan model installation details are 
provided subsequently. 


Code Output Results 

Figures 3 - 5 present a sample of the WICS rode 
output from the NTF for a seinispan model. Figure 3 
gives a row-wise com]>arison of the tared signature 
with the WICS-fit wall signature for a single test point 
at M = 0.2,n„. m = 24° for a semispan model. The 
symbols represent measured, tared wall perturbation 
velocities and the solid line represents the global least 
squares prediction of the code. (One* measure of va- 
lidity is how well these two match.) Mean primary 
corrections arc* displayed in Figure 4. These correc- 
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Fig. 3 Wall Signature of a Low Aspect Ratio-High 
Lift-Semispan Model at o„ nr = 24° 

General Implementation Requirements 

To implement the WICS code, 4 * several tasks must 
first he completed. The first task is to ensure the facil- 
ity has an adequate? number of static pressure? orifices 
on the wall of the test section. (The 1 meaning of ad- 
equate will be discussed later in this paper.) Next, a 
perturbation velocity database 1 (PVD) must be gener- 
ated using a panel method or the method of images. 
This PYD is used in the matching of the tared wall 
signature to determine singularity strengths, and to 
generate the field of interference velocities so that cor- 
rections can be computed. 

Once the orifices have been installed, an “empty t mi- 
ner' calibration must be performed. This calibration 
is a function of independent test parameters ( e.g. y to- 
tal pressure and Mach number) and support system 
attitude (e.g., pitch and sideslip angles). For semis- 
pan models, the empty tunnel calibration is the tunnel 
geometry minus the model. For fullspan models the 
model support system is included, and the sting or 
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b) Incidence Correction 

Fig. 4 Mean Primary Corrections for a Low As- 
pect Ratio-High Lift-Semispan Model 

tions are shown for an entire pitch polar (or group of 
data points). Blockage corrections (c) are averaged in- 
terference velocities along the fuselage centerline, and 
upwash (incidence) corrections (a, ) are averaged along 
the wing three-quarter-chord. Corrections to Mach 
number (AM) and dynamic pressure (AQ) are de- 
rived from the blockage factor. Coefficient corrections 
for the entire polar are plotted in Figure 5. Note for 
the test point of a — 24° in Figure 3 that a blockage 
induced correction of AM = 0.0025 and a significant 
incidence correction of o, = 1.33 are determined as 
shown in Figure 4. Correspondingly large corrections 


•I or 12 


American Institute of Aeronautics and Astronautics Paper 2001 2421 



0 008 



0 006 
0 004 
0 002 
0 000 
-0 002 
-0 004 
-0 006 
-0 008 
•0010 
.-0 012 


; 0.0010 
5 0 0008 
\ 0 0006 
| 0 0004 

j 0 0002 
j 0 0000 

; -o 0002 

; -0 0004 
: -0 0006 
: -0 0008 
' -oooio 


Fig. 5 Coefficient Corrections for a Low Aspect 
Ratio-High Lift-Semispan Model 

for the coefficients an* also observed (Figure 5). Since* 
the WICS code calculates the* interference corrections 
for specified reference planes in a volume surrounding 
the model, contours of tin* correction variation can he 
plotted. 


WICS code responds to the "soft wall" as opposed to 
the "hard wall" assumptions made in classical theory. 
If there is no correction for wall effects, a pseudo- 
Reynolds number effect may he embedded in the wind 
tunnel data. 

Mach Number 

According to classical wall interference theory, 
the blockage 1 factor c should scale* according to the 
Prandtl-Glauert compressibility factor H. Figure* 7 
shows the blockage factor e>ve*r a data run for several 
Mach numbers. It can be seen that tin* variation e>f .57 
at constant lift coefficient is not a linear function of 
Mach number; in fact., it is a function of order greater 
than or equal to three. Insufficient data e*xist to further 
evaluate this phenomena: therefore, this is an area for 
further study into tlx* tunnel How field to determine 
the source of this phemonena. 


Aerodynamic Parameter Sensitivities 

The purpose of this section is to understand how 
the WICS code responds to the independent variation 
of Reynolds Number and Mach Number. (Note that 
NTF also allows the independent variation of dynamic 
pressure, but no data of this type were available for 
use in this study.) 

Reynolds Number 

As previously stated, boundary pressure methods 
allow inclusion of actual tunnel response in the com- 
putation of wall effects. This is due to a response to 
the pressure change due to the boundary layer height 
changing, which is highly dependent on Reynolds num- 
ber. This can easily be seen by plotting the total 
blockage factor, c 7 . against the uncorrected lift co- 
efficient., Ci Jtinv . for varying Reynolds numbers at a 
Mach number of 0.2 (set* Figure G). Blockage lev- 
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Fig. 6 Total Blockage Variation with Reynolds 
Number (M — 0.2) 

els arc* inversely proportional to Reynolds number, 
due to higher Reynolds numbers causing thinner wall 
and model boundary layers, corresponding to an effec- 
tive* change in tunnel cross-sectional area. Thus the 



Fig. 7 Total Blockage Variation with Mach Num- 
ber 

Support System Sensitivities 

The purpose* of this section is to discuss the* potential 
error caused by support, systems. Model supports for 
both fullspan and semispan models will be discussed. 

Fullspan Support Systems 

A major concern with the use* of the WICS code* 
in tunnels is the cost of the* fullspan support system 
calibration. While* this is necessary to obtain modeT 
only wall interference* corrections, it can be e*xtremely 
expensive* and on occasion unfe*asible*. I lbrich 1 1 has 
developed an alternative method which doc*s not re- 
quire a support syste*m calibration for sting mounted 
models. His alternative* method re*places the support 
system solid blockage source singularities by chains 
of point doublets. These chains an* placed along the* 
support system centerline and are* weighted aceoreling 
to the cross-sectional area distribution. This method 
does not currently work with post mounted models. 

In terms of the code* sensitivities, if the* support sys- 
tem calibration can be* shown to be n<*gligible. it could 
be* ignored and the true* empty tunnel data could be* 
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used in its place. Unfortunately, the interference in- 
duced by almost all fullspan support systems will cause 
a noticeable error in the corrections if they are ne- 
glected. This error in the corrections was also noted 
by Murthy. 15 

Semispan Standoff 

Questions have been raised concerning the calcula- 
tion of wall interference corrections with a semispan 
model standoff support system present. The standoff 
contributes to the blockage and lift interference in the 
tunnel but is non-metric with respect to tin* balance. 
Since the WIGS code incorporates data from both the 
wall pressures and the balance, an inherent, inconsis- 
tency in the method of correction is introduced. 

For a fullspan model the “empty tunnel’’ wall signa- 
ture includes the support hardware, unlike the semis- 
pan technique. The pressure signature of the standoff 
cannot, easily he identified and removed from the wall 
signature, if at. all. The flow field produced bv test- 
ing only the standoff is radically different than with 
the model installed creating an inconsistency dilemma. 
Two solution approaches are presented. First, depend- 
ing on the type of testing, the standoff blockage may 
be neglected. Second, the lift on the standoff can be 
calculated so the bookkeeping of interference becomes 
consistent. 

Studies have been performed to simulate the lift in- 
terference' of the standoff. A transport semispan model 
was tested in the NTF. Centerline fuselage pressures 
were measured. These pressure's where integrated in 
the chord-wise 1 direction around the fuselage to obtain 
a twodimensiona] lift coefficient for each angle of at- 
tack tested. Milholen and Chokani 16 have shown that 
the chord Cp distribution can be assumed to be ap- 
proximately constant from the fuselage' centerline to 
the wall (Gc., standoff root to tip) if the standoff is a 
two-dimensional extension of the fuselage centerline. 1 ' 
These assumptions allow the span-wise integration of 
the chord lift across the width of the standoff, thus 
yielding a three-dimensional lift coefficient. 

CL.„=jc P d{L)4jfc) (5) 

Here ~ is the percent chord and is the height of the 
standoff normalized by the semispan. Figure 8 shows 
the magnitude of the calculated lift increment due to 
the standoff. 

Test data were then corrected with and without the 
standoff lift increment. Figure 9 shows this correc- 
tion comparison. Not including the lift of the stand- 
off causes an over prediction in the blockage. This 
happens because the value of lift is changed and the 
corresponding singularities also change. The result of 
this increase in lift is that more lift interference is sub- 
tracted from the wall signature, thereby, causing the 
bloc kage interference to be lower. Increase's are seen in 
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Fig. 8 Integrated Three-Dimensional Lift of a 
Two-Dimensional Semispan Transport Standoff 

the incidence correction due primarily to the increase 
in lift. 

Walker 18 and Walker et a/. 19 have determined the 
sensitivity of the code* to error in the lift measure- 
ment. Since the method of correction for the presence 
of the semispan standoff entails only a incremental 
adjustment to the measured lift, this can be treated 
as causing an error to occur in the lift measurement. 
Given this logic, correction for the presence of the 
semispan standoff may be safely neglected when the 
lift generated by the standoff falls within the general- 
ized correction accuracy requirements given in Table 1. 



Fig. 9 Effect of Standoff Lift on Primary Correc- 
tions 

Other Sensitivities 

In an effort to reduce the influence of outlying data 
on the global least squares solution, a data rejection al- 
gorithm was implemented in the WIGS code. This sec- 
tion addresses the standard data rejection algorithm of 
the WIGS code. The primary topic of discussion will 
be the algorithm’s interaction with the sensitivities. 
Two different types of improvements are suggested. 
The first involves a more statistically meaningful rejec- 
tion algorithm, and the second involves a more robust 
matrix solver that is not easily biased by outlying data. 

Standard Data Rejection Algorithm 

Currently, the algorithm rejects data based on two 
criteria. The first test performed is an absolute range? 
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chock. If wall signature data exceeds 20% of the cor- 
responding empty tunnel calibration, those data are 
not used in the correction calculation. The second 
test attempts to locate and remove' outliers. It does 
so by calculating the average difference between the 
tared wall signature and the wics fit. Data whose ab- 
solute' difference is greater than 3 times the RMS of 
the average' difference are removed from the* calcula- 
tion. Essentially, the* classical form of the standard 
deviation is used and data which lie' more than three 1 
standard deviations from the mean are' rejected. The 
sample variance 1 is computed as follows: 

X! ~ 11 fit) 

- x-r-1 m 

where' is the measured in tunnel perturbation 

velocity. Uf, f is the calculated least squares fit of u tlini , 
N is the 1 number of active 1 ports, and P is the number 
of j)aramet('rs or dimensions in the solution process. 
Note* that in the original version of WICS the value 
of P was set to zero, but this had very little 1 effect 
on the calculated variance due* to the* large value of 
A (2G0-290 in NTF). Once 1 these 1 tests are 1 performed, 
the calculation procedure* is only allowed to continue 
if more 1 than GO port values remain in the 1 calculation. 
The 1 flow field is them split into three* parts: upstream 
of the* inoek'l, around the model, and downstream of 
the model. The active wall pressure port distribution 
is calculated and reported as a percentage of the total 
active ports in each area of the 1 flow field. 

The data rejection algorithm was evaluat.e*d follow- 
ing some 1 concern that valid data were b(*ing discarded. 
Normally this would not be an issue* sinc e* the 1 WICS 
code 1 computes a global least squares fit, erf several hun- 
dre'd wall pressures using the method of singular value* 
decomposition to determine the blockage* singularity 
strength and thereby the blockage. So the* loss of a 
relatively small numbe*r of valid data points basically 
has little to no e'ffect. on the solution. However, the 
port placement sensitivities were defined by Walker 18 
and \Valke*r ct a/.. 19 and regions of critical ni(»asure»- 
ments were 1 established. Ports around the* model and 
wake* singularit,i(*s are implicitly weighted more 1 heav- 
ily than those* furthe*r upstream or downstream. The 
implicit weighting approximately varies as the inverse* 
of the* square of the distance* from the singularity to 
the* wall port for sources (i.e., influence* decays quickly 
with distance 1 ). The 1 implicit weighting for doublets 
decays as the inverse of the cube of the radial distance 
from the singularity to the 1 wall port. Accordingly, 
the loss of valid measurements in regions of relatively 
strong influence can adversely affec t the sensitivities 
of the output, corrections. 

The ('fleet of the data rejection algorithm on the 
wall signature is shown in Figure 10. This figure* con- 
tains four representative rows of data near the mode 1 ! 
(Figure’s lUa-lOd) and a break out of wall signature 


components for a row (Figure 10e). The filled circle's 
represent data that were* retained afte»r execution of 
the* data rejection algorithm. Data that were rejected 
are 1 indicated by the* open curie's. The* solid line is the 
WICS calculated fit of the* retained data. The* code* 
was again executed but forced to keep all wall data. 
This yielded a fit of all the data or the* total fit, which 
is denoted by the dashed line. 
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Fig. 10 Data Rejection Algorithm Effect on the 
Wall Signature of a Low Aspect Ratio-High Lift- 
Seniispan Model at = 24° 


As can be* seen from Figure lOe, resolution of the 
wall signature 1 is lost when wall data arc* lejerted. This 
is evidenced by the* straight line* connection among the 
rejected data in the retained data curves. Since the 
only change between the two case's is the inclusion of 
wall data that would have otherwise been rejected, the 
doublet, strengths representing the lift component of 
the wall signature* remain unchanged. Even though the 
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strengths remain the same, the lift signature changes 
when all the wall measurement, points are included 
because the wall flow field has been resolved in finer 
detail. This in turn has a direct impact on the 1 least 
squares solution of the solid and wake blockage singu- 
larities since a different value 1 for the lift interference 1 
has been subtracted at the previously rejected ports. 
The 1 effect of the data rejection algorithm on the pri- 
mary corrections is shown in Figure 11. The overall 
result is that rejection of the valiel data near the model 
singularities caused the code to over predict blockage 
and. in turn, the correction te> model incidence. The 
conclusion is drawn that when valid data in legions 
of high sensitivity are discarded a true global solution 
cannot be determined. 
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Fig. 11 Data Rejection Algorithm Effect on the 
Mean Primary Corrections of a Low Aspect Ratio- 
High Lift-Semispan Model 


Improved Data Rejection Algorithm 

Taking into account the regions of high sensitivity 
near the model singularities and the W IC S mathemat- 
ical model, a modification to the calculation of the 
standard deviation has been developed to improve the 
data rejection algorithm. The WICS mathematical 
model contains two desc riptive variables with no con- 
stant term 



where the elements of A are determined from the 
PVD, the elements of B are determined by subtracting 
all but the blockage component from the measured wall 
data, and the elements <7* and rr** arc 1 the solid vol- 
ume and viscous wake blockage singularity strengths, 
respectively. 

One of the major issue's is that the residuals — 

Ufa) are not normally distributed. From this, it can 
be inferred that ujn is not the true mean of the wall 
data. This fact is quite possibly due to the low order 
of the mathematical model of the flow field (i.e., there 
arc 1 insufficient degrees of freedom in the mathemati- 
cal model to allow for the solution of the true mean, 
even if no error is present in the data). Given this 
reasoning, it is easily understood why the data in Fig- 


ure 10 were rejected. Use of ujn as the sample mean 
has biased the calculation of the standard deviation, 
t.herelw under-inflating the confidence interval set to 
reject data. Now that, it is understood that. uf i( is a 
biased estimator for the sample 1 mean of u timi . a better 
estimate of the confidence interval can be constructed 
by assuming that Uf,t is a low-order estimate ol the 
sample mean. This means that tin* corresponding de- 
grees of freedom in the 1 solution should be 1 reduced such 
that a proper inflation of the confidence interval can 
be achieved. The reduction in the number of degrees 
of freedom can be achieved by considering the implicit 
weighting of the singularity strengths based on the* in- 
dial distance 1 from the model singularities to the 1 wall 
pressure 1 ports. 

Define r* as the minimum radial distance 1 from a 
singularity to a wall pressure port. Then a weighting 
(IF) of the influence each wall pressure port has on the 
solution can he determined for a source by 



and for a doublet by 

r :i 

W* = ^ (9) 

r i 

The number of effective wall orifices (A«. ) in the tun- 
nel is calculated by summing the influential weights of 
each orifice 

A’„ ='£ W ' (lt)) 

> = l 

Philosophically, W should be determined for each sin- 
gularity solved for, g * and o'**, and the nuinbei of 
effective wall ports, N w is the difference between the 
union and intersection of the sets. 


n w = w itrm u w inmt - n W... ( n ) 

where \\\ m and W inmt denote the set of weights for 
each wall port based on the distance to the respective 
singularity. The sample variance is then calculated by 

•> _ T. ( U tuuj — ( 12 ) 

N w -P-1 

where P is the number of parameters or descriptive 
variables in the mathematical model, for the W ICS 
code P = 2. Use of this method allows a more 
proper inflation of tin* confidence interval such that 
data which are truly outlying and have a potential of 
biasing the solution will have a much better probabil- 
ity of detection and removal. 

Solutions using Robust Regression Techniques 

As an alternate solution to the problem of outlying 
data, robust regression methods have been evaluated 
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method is given by 


for use in solving the solid volume and viscous wake 
blockage singularity strengths. The methods of Least 
Median of Squares (LMS) and Least Trimmed Squares 
introduced by Rousseeuw and Leroy 20 were analyzed 
as potential matrix solvers for the WTCS rode. Use 
of these techniques allows the solution some level of 
insensitivity to error in the measured data. Then 1 are. 
however, some disadvantages to using these techniques 
with this type of modeling or simulation. These will he 
discussed after a presentation of the basic methods. To 
facilitate an easier comparison of the regression tech- 
niques, a brief description of the least squares method 
is included. 

Consider a general linear model with P explanatory 
variables, one response variable, and N data points. 
This model can be expressed as 

Y\ x , = D V)(r X/. xl 4- E v > i (13) 

when 1 Y is the response variable, D is the matrix 
of explanatory variables. X is the solution vector of 
coefficients, and E is the vector of errors, which is clas- 
sically assumed to be normally distributed with mean 
of zero and unknown variance* o ~ . The relation for t he 
solved estimation is expressed as 

Y\ *i = D iVxJ . X/» xl (14) 

where the hat denotes an estimated parameter or vari- 
able. From these expressions the* matrix of residuals is 
defined as 

HjYxl — Y Vx j — Y a x | (15) 

The method of least squares (LS) then corresponds to 

jY 

Minimize fij based on X (16) 

j= i 

An important topic of discussion is the notion of a 
breakdown point. Basically, the breakdown point of 
a regression method corresponds to the percentage of 
contamination allowed in the data set such that if an 
exact fit exists the method is able to resolve it. For 
the LS method, its breakdown point is 0 % . In other 
words, if there* exists any data contamination the LS 
method will not resolve an exact fit. 

Unlike the mean of a set of data that is influenced 
easily by any outlying data, the median is very robust. 
This is tlu* basis for the LMS method which can be 
expressed as 

Minimize* median,/? 2 based on X (17) 

The breakdown point of the LMS method is 50 %., 
which is the highest possible value*: unfortunately, 
Rousseeuw has shown that LMS has an abnormally 
slow convergence rate*. In an effort to increase the con- 
vergence* rate, he* introduced the LTS method. This 


h 

Minimize (/?j) v based on X (18) 

where (/? 2 ) |V < < (/? 2 ) vv are ordered squared 

residuals. The LTS me*thod is similar to LS with the* 
exception that the largest residuals are* not. included in 
the calculation. This is the so called trimming of tin* 
LS method. The breakdown point of the* LTS method 
depends on the value* chosen for h. A breakdown point 
of 50% is achieved when /; is approximately n/2. 

The* matrix solver in the WTCS code was replaced 
with both a LMS and LTS solver for comparison. The* 
algorithms were* adapted from tin* Program for Robust 
Regression (PROGRESS) described by Rousseeuw and 
Lerov J0 and improved by Rousseeuw and Hubert. 2J 
The breakdown points for both methods we*re* set to 
25 ( /f for reasons discussed below. The computation 
lime of the* full LMS and LTS methods is extremely 
slow compared to the WTCS LS calculation. Since 
speed of computation was one* of the* main reasons for 
the choice of the* W ICS code*, a s(*ve*re increase in com- 
putation time is not desirable*. Fortunately. Rousseeuw 
and \an Driessen - ~ have* developed a LTS approxi- 
mation algorithm (FAST-LTS) that executes at ap- 
proximatc'ly th<> same speed as the standard WIC’S LS 
method. 

Figure* 12 details a comparison of the wall signature 
computed from the various methods, and Figure* 13 
shows the insulting primary corrections. At this time* 
no final conclusions have* been drawn concerning the* 
use of robust regression methods in the* WTCS code. 
There* were, however, some* important observations. 
Since the* mathematical mode*! in the WTCS code (refer 
to Equation 7) lias been reduced to a simplistic re*pr<»- 
sentation of the* flow field, use* of high-breakdown point 
regression estimators may also result in the insensitiv- 
ity to certain local physical fluid dynamic phenomena 
due the large residuals produced when the mathemat- 
ical model cannot be forced to fit the measured data. 
The combination of this insensitivity effect and the 
lack of resolution of the wall signature* far downstream 
of the model and wake* in NTF resulted in the selec- 
tion of a 25% breakdown point for the* LMS. LTS, and 
FAST-LTS methods. The high-breakdown point r<*- 
gression methods an* anticipated to perform better for 
the WTCS implementation at the NASA Langley 14-by 
22-Foot. Subsonic Tunnel (14 x 22) due* to a proposed 
increase in the resolution of the* wall signature far down 
stream of the* model. Implemientatiem of the WTCS 
code* at. 14 x 22 is discussed by Iyer and Everhart. 2:i 
Another important observation was the* strong agree- 
ment between all the methods in regions where little* 
or no separated flow occurred. This is indicative of 
the greater uncertainty involved in the computation of 
flow fields with regions of separated flow. 
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Conclusions 

This paper makes several general conclusions con- 
cerning the WICS code. The code appropriately re- 
sponds to change's in Reynolds number. Some ques- 
tion arose as to the code’s response to Mach iiumbei 
changes at constant Reynolds number, but additional 
data are required to further investigate the true re- 
sponse. 

Fullspan support systems almost always create 
enough blockage that they must be accounted for by 
some method. For semispan models using a stand- 
off, a correction to the lift based on the standoff is 
necessary if the standoff lift is predicted to generate 
more lift than the generalized correction accuracy re- 
quirements allow. Typically this entails a standoff lift 
greater than one-haif of the absolute correction accu- 
racy requirements or 0.2% of the model generated lift. 

It was also concluded that there were circumstances 
in which the standard data rejection algorithm of the 
WICS code rejected valid data which had an influ- 
ence on the sensitivities of the corrections. To avoid 
this problem an alternate rejection algorithm was pro- 
posed. In an attempt to completely circumvent the 
rejection algorithm an evaluation of robust regression 
estimators was begun. No final conclusions were made 
for the use of a robust regression method in the WICS 
code, hut several findings were presented. More anal- 
ysis should he performed to properly assess the use of 
these methods in this type of engineering data analy- 
sis. 
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Fig. 12 Comparison of Solution Algorithm Effects on the Wall Signature of a Low Aspect Ratio-High 
Lift-Semispan Model at a = 24° 
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Fig. 13 Comparison of Solution Algorithm Effects on the Primary Corrections of a Low Aspect Ratio- 
High Lift- Semispan Model 
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